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ABSTRACT 

Context. Theoretical arguments suggest that dust grains should grow in the dense cold parts of molecular clouds. 
Evidence of larger grains has so far been gathered in near/mid infrared extinction and millimeter observations. 
Interpreting the data is, however, aggravated by the complex interplay of density and dust properties (as well as 
temperature for thermal emission). 

Aims. Direct evidence of larger particles can be derived from scattered mid-infrared (MIR) radiation from a molec- 
ular cloud observed in a spectral range where little or no emission from polycyclic aromatic hydrocarbons (PAHs) is 
expected. 

Methods. We present new Spitzer data of L183 in bands that are sensitive and insensitive to PAHs. The visual extinc- 
tion Ay map derived in a former paper was fitted by a series of 3D Gaussian distributions. For different dust models, 
we calculate the scattered MIR radiation images of structures that agree with the Av map and compare them to the 
Spitzer data. 

Results. The Spitzer data of L183 show emission in the 3.6 and 4.5 /im bands, while the 5.8 fira band shows slight 
absorption. The emission layer of stochastically heated particles should coincide with the layer of strongest scattering of 
optical interstellar radiation, which is seen as an outer surface on I band images different from the emission region seen 
in the Spitzer images. Moreover, PAH emission is expected to strongly increase from 4.5 to 5.8 /im, which is not seen. 
Hence, we interpret this emission to be MIR cloudshine. Scattered light modeling when assuming interstellar medium 
dust grains without growth does not reproduce flux measurable by Spitzer. In contrast, models with grains growing 
with density yield images with a flux and pattern comparable to the Spitzer images in the bands 3.6, 4.5, and 8.0 ^m. 
Conclusions. There is direct evidence of dust grain growth in the inner part of L183 from the scattered light MIR 
images seen by Spitzer. 

Key words. ISM: dust, extinction - ISM: clouds - ISM: individual object : L183 - Infrared: ISM - Radiative Transfer 
- Scattering 



1. Introduction 

Cold and dense cores of molecular clouds are commonly 
considered to be a favorable place for the grow th of dust 
particles (see, e.g.. lOssenkopf fc Henninglll994l ). Evidence 
that grains are larger than commonly observed interstel- 
lar medium (ISM) grains have been gathered in different 
wavelength regions. The thermal emission of dust grains 
in the far-infrared (FIR) and millimeter (mm) wavelength 
range shows a change in the index of the spectral den- 
sity distribution. Based on simple assumptions about the 
dust opacity law and the spatial and thermal structure, 
this is interpreted a s indirect evidence of larger grains (e.g. 
Stepnik etall 120031; iBeuther et all 12004; iRidderstad et all 
2006t iKiss et al.l 120061 : ISchnee et all I2008D . In contrast, 
Nutter et al. ( 20081 ) can model the Taurus Molecular Ring- 
without any assumption about grown particles. In general, 
disentangling spatial variations in dust properties, density, 
and temperature is a challenge given the many parameters 
needed to describe the often clumpy or filamentary cloud 



structures. This is why, in view of this ambiguity, we call 
this indirect evidence. 



[Chap man et al] (120091) . iMcClurd (I2009I) . and 



ICha^ man fc Mundvl ((2009) have presented mid-infrared 
(MIR) extinction curves based on Spitzer and ground- 
based JHK band observation of stars behind molecular 
clouds. They find that the extinctio n curve c hanges 
as a function of increasing extinction. iMcClur el (|2009( ) 
investigates the ice features and concludes that a process 
involving ice is re sponsible for th e chang ing shape of the 
extinction curve. IChapman et all (|2009l ) find that, with 
increasing extinction, the extinction law is better fitted 
by models incl uding larger grains. S i miliarl y, for the cores 
investigated in IChapman fc Mundvl (|2009l ). the averaged 
extinction law from 3.6 to 8 /im is consistent with a dust 
model that includes larger grains. 

Extinction measurements at optical to MIR wavelengths 
have the advantage that the temperature of the dust grains 
no longer influences the results. The disadvantages are (i) 
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the method requires a radiation source like scattered light, 
diffuse emission from small grains or molecules, or back- 
ground stellar or cxtragalactic emission; (ii) the dust ex- 
tinction properties enter exponentially in the observed ex- 
tinction of the external radiation so that the source flux, 
if stellar or extragalactic, must be strong enough to be de- 
tectable. 

Aside from the thermal emission of the grains and their 
absorption, they also efficiently scatter radiation. The ex- 
ploration of the physical properties and spatial distribution 
of dust particles via the observation of their scattered light 
was initiated long ago. Referen ces to early paper s abou t 
this topic can be found, e.g., in ISandell fc Mattilal (|l975h . 
Common ISM grains with sizes around 0.1 fim efficiently 
scatter optical and near-infrared (NIR) radiation so that 
images observed at these wavelengths can be used to study 
the spatial structure of the dust, as well as its properties. 
This method has widely been used in star formation re- 
search, e.g. to analyze the spatial 3D structure of circum - 
stellar disk candidate s (see, e.g., ISteinacker et al.1 [2006( 1. 
iNakajima et a l. (2003) interpret the surface brightness seen 
in J, H, and Ks ima ges of the Lupus 3 da r k clou d as starlight 
scattered by dust. iFoster fc Goodman! (|2006h propose to 
use this "cloudshine" as a general tracer of the column den- 
sity of dense clouds and presented NIR scatt ered light im- 
ages o f the Perseus molecular cloud complex. iPadoan et al.1 
(|2006h present a corresponding method and test it by calcu- 
lating NIR scattered light images from the results of magne- 
tohydrodynamic turbulence simulations and by comparing 
the resulting column density distributions with the origi- 
nal simulation data. They argue that the method can be 
applied to filamentary cloud structures being illuminated 
by a homogenous radiation field in the range of visual ex- 
tinction between 1 and 20 mag. The NIR cloud shine of 
a Coro na Australis cloud fi lament was investigated in this 
way bv lJuvela et~aTl (|2008l ). 

It seems promising to observe scattered light from cores 
at longer wavelengths, because the lower extinction allows 
the inner parts to be traced. However, the scattering cross 
section of the ISM dust particles drops with wavelength 
once the wavelength is much larger than the particle size. 
Scattered light from cores at MIR wavelengths therefore 
has not been reported so far and has been overlooked un- 
til now despite submillimcter observations often being ex- 
plained by grain growth to sizes for which MIR scattering 
may become visible. Additionally, polycyclic aromatic hy- 
drocarbon (PAH) emission starts to dominate the emission 
pattern, at most MIR wavelengths, preventing the identifi- 
cation of scattered light. 

It turns out that indeed the scattered light image anal- 
ysis at MIR wavelengths opens up a window to trace bigger 
grains: (i) they scatter radiation efficiently at wavelengths 
for which the smaller grains already enter the regime of 
Rayleigh scattering where absorption dominates scattering; 
(ii) the optical depth is reduced compared to the depth at 
NIR wavelengths so that denser regions can be traced with 
a better chance of seeing grown particles. 

Observations with the Spitzer telescope allow us to trace 
the relevant wavelengths and to check the existence of pos- 
sible PAH emission by using a combination of filters, some 
including strong emission features and some excluding all. 
Nearby starless low- to mean-mass clouds are potential can- 
didates for observing such an MIR cloud shine. 



L183 (also known as L134N) is such a starless dark cloud 
hig h above the G alactic plane (36°), hence close to us (110 
- Francdll989h. Its mass has been estimated to ^80 Mq 



CEa gani et al 12004 hereafter Paper I). It contains an elon- 
gated ridge of dense material with two prestellar cores and 
the main prestellar core has a peak extinction of Ay ~ 150 
mag with a dust temperature close to 7 K (Paper I) . It was 
subsequently sho wn that the gas is thermalizcd with the 
dust in the core (|Paeani et al.ll2007D . 

The Spitzer observations of L183 are summarized in 
Sect. [H In Sect. [31 we fit the Ay map presented in Paper 
I with a series of 3D density basis functions. The results of 
a scattered light image modeling are presented in Sect. [4] 
based on a dust growth model and compared to the differ- 
ent Spitzer images of L183. In Sect. (5[ we summarize our 
findings. 

2. Spitzer observations of L183 

To describe different parts of the cloud, we refer to the Ay- 
map presented in Paper I (see their Fig. 1 for a detailed 
comparison), although it is a projected view on the spatial 
cloud structure. Throughout this paper, we use the term 
outer cloud region for the parts of the LI 83 cloud where 
Ay< 15, (banana-shaped) central part for 15 < Ay< 90, 
and central core for the densest region with Ay> 90, re- 
spectively. 

The Spitzer IRAC (3.6 to 8 ^m. iFazio et~aIll2004D ob- 
servations of LI 83 were obtained as part of the Cycle 1 
GTO program 94 (PI: Charles Lawrence). The observations 
were broken into two epochs, the astronomical observa- 
tion requests (AORs) are both centered at 15h54ml5.00s - 
2d55m0.0s, they were obtained on 23 and 25 of August 
2005, and the AORKEYs are 4921600 and 4921856 in 
chronological order. Each AOR was constructed with the 
same strategy: at each map step, three dither positions were 
observed, each with high-dynamic-range exposures of 1.2 
and 30 seconds. 

To build the full mosaics of L183, we started with the 
Spitzer Science Center (SSC) pipeline-produced, basic cal- 
ibrated data (BCDs), version S14.0. We ran the IRAC 
artifact mitigation code written by S. Carey and avail- 
able on the SSC website. We constructed a mosaic for 
each epoch from the corrected BCDs using the SSC mo- 
saicking and point-source extraction (MOPEX) software 
([Makovoz fc Marleaul[2005[ l. with a pixel size of 1.22" px" 1 , 
very close to the native pixel scale. 

The inner parts of 3.6, 4.5, and 8 /im mosaics are shown 
in Fig. [TJ One can clearly see a banana-shaped structure 
(see also Sect. 14. 2[ ) visible in emission at 3.6 and 4.5 /im and 
in absorption at 8 fim. The spatial structure resembles well 
the shape of the densest region of the cloud complex visible 
in the Ay map presented in Paper I. The 5.8 /im image 
shows horizontal stripes and a strong vertical gradient in 
signal-to-noise ratio (S/N) so that only horizontal pieces of 
the upper part of the image contain reliable information on 
the flux. For this reason, we do not show the image here, and 
do not include it in the image modeling. The upper part of 
the central banana-shaped structures is seen in absorption 
in the reliable stripes. 

To show the magnitude of the emission and absorp- 
tion features, we present in Fig. [2] cuts of the IRAC maps 
through the central core (declination -2°52'48") seen in ab- 
sorption on all images. As the cut lies in a stripe of the 5.8 
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Fig. 1. IRAC 3.6, 4.5, and Sfim long exposure mosaics. Superposed on each image are the Ay — 5 & 10 mag contours 
from the Ay map of Paper I. The white crosses in the 8 image mark the region that has been modeled in this paper. 
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Fig. 2. Cuts from the IRAC maps (Fig. [T]) at constant dec- 
lination (-2°52'48") through the dust peak in all 4 bands. 
The dotted line indicates the position of the dust peak. 
Excess emission is seen at 3.6 and 4.5 /zm while absorption 
dominates for the 5.8 and 8.0 /xm ranges. 



(im image with an acceptable S/N, we also include the flux 
in this band. The gray thin lines indicate an averaged exter- 
nal flux. From the cuts, the increase by almost two orders 
of magnitude in the external flux is visible when changing 
from 3.6 to 8 /im. 

The thermal emission of dust particles which are large 
enough to be considered as black body emitters cannot pro- 
duce a measurable flux gi ven the low temperature of a few 
ten K (jPagani et al.ll2004h . Even if all dust particles in L183 
would have a temperature of 40 K, the specific flux density 
measured on Earth at A = 8 ^m would be < 7 x 10 -29 Jy/sr 
(assuming 0.1 ^tm-sized grains and a gas-to-dust mass ratio 
of 100), and lower for the shorter IRAC band wavelengths. 

The emission in the Spitzer wavelength range could also 
be caused by smaller particles that are stochastically heated 
like PAHs. Figure [3] illustrates the presence of PAH emis- 
sion features in IRAC bands by showing the instrumen- 



tal response of IRAC (gray areas) for the different filters 
as a function of wavelength. The transmission of radiation 
thr ough a layer o f dust is plotted as a black dashed line 
fsee lDrain e 2003a. b]Q). The light-gray (red in the electronic 
version) line shows the expected location of PAH features 
(iDraine fc Lill2007t ). 

Explaining the observed emission in the 3.6 and 4.5 
/im bands as the emission of stochastically heated parti- 
cles (PAHs or small dust grains) is difficult for two rea- 
sons. First, PAHs and small grains need an exciting UV to 
optical radiation field. As dust particles have the same or- 
der of magnitude in absorption and scattering cross section 
at these wavelengths, the scattered light images at optical 
wavelength give a hint of the morphology of the layer where 
most of the outer exciting UV t o optical radiation is ab- 
sorbed. As seen, e.g, in Fig. 4d of lJuvela et all (|2002t) . this 
layer resembles the shape of the Ay=l-layer in the Ay map 
closely (Paper I). In contrast, the emission pattern resem- 
bles the shape of the inner banana-like structure while no 
emission can be seen with the shape of the outer layer. 
Even if holes would make it possible for exciting radiation 
to reach inner parts, a substantial part of the PAH emission 
should come from the outer layer which is not what is seen. 

For completeness, we also investigated the possibility 
of holes. We performed detailed 3D ray-tracing for a large 
sample of structures that are consistent with the column 
density information contained in the Ay map. The results 
are presented in Appendix [Al Based on 3D spatial structure 
variations and the analysis of r-iso-surfaccs, we conclude 
that only a small fraction of all spatial structure models 
would allow for PAH emission near the central absorption 
peak. It therefore seems unlikely that the cloud complex has 
one of these rare configurations with a hole deep enough to 
let the interstellar radiation field (ISRF) reach the central 
part so that PAHs could be excited to produce the observed 
inner emission. Moreover, the illumination through a single 
hole would not produce a smooth, even illumination of the 
entire inner structure as seen in the 3.6 /*m image. 
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Fig. 3. The gray areas are the instrumental response of 
each IRAC filters as a function of the wavelength. The 
light gray (red in the electronic versio n) line is a model o f 
infrared emission from dust grains bv lDraine & ~L1 (120071 ) 
(q_PAH=4.58% U=l). T he black da shed lin e is a model of 
the dust extinction from iDrain 1 (see l2003allbl lcl). The mod- 
els are unitless, and have been scaled to fit in the figure. 



Second, based on an integration of the gray/red curve 
in Fig. [51 the PAH emission in the 5.8 fim band should be a 
factor of about 10 stronger than the emission in the 4.5 (im 
band. The improved dust transmission in the 5.8 /an band 
will amplify this effect. In turn, the external radiation at 5.8 
li m is about a factor o f 10 stronger than at 4.5 ^m. Based 
on iFlaeev et all (|2006l) . we estimated the flux one would 
expect if the excess flux at 4.5 ^m was due to PAHs. We 
found that, in all reasonable cases, the 5.8 /im flux should be 
seen in emission (about 1 to 2 orders of magnitude stronger 
than at 4.5 /im) and not in absorption as presently observed 
(see cuts in Fig. [5]). 

These two arguments suggest that emission from 
stochastically heated grains likely fails as an explanation 
of the emission structure. We therefore explore in this pa- 
per whether the emission seen in the 3.6 and 4.5 /xm filters 
can be interpreted as stellar light scattered by the dust par- 
ticles in the cloud gas, and what dust properties would be 
needed to qualitatively reproduce the images. 

3. Ay map modeling 

Extinction maps average the 3D density distribution of the 
molecular cloud over the line of sight (LoS). This adds 
an important source of ambiguity to the interpretation of 
the extinction map: each observed extinction value in the 
plane of sky agrees with a whole range of density distri- 
butions and dust properties. As a result, another source of 
information is needed to remove this ambiguity. This can 
be, e.g., a measurement at FIR/mm wavelengths. Multi- 
wavelength modeling of images from molecular clouds and 




Fig. 4. Inner 66000 AU of the A v map from Paper I (left) 
and the fitted map using 3D Gaussian clumps (right). 

cloud cores incorporating the hidden third dimension has 
not been used widely so far , but i s possible in the mean time 
(see, e.g.. ISteinacker et al.ll2005l ). Another source of infor- 
mation are scattered light images. They depend on the 3D 
spatial structure of the cloud, on the dust properties, and 
the external 3D radiation field. 

In this work, we combine the information about L183 
from the extinction map with the scattered light images ob- 
tained by Spitzer. To obtain a 3D dust density data cube 
required for the scattering calculations, we fit the available 
Ay map of L183 (Paper I) with a series of density basis 
functions flexible enough to deal with the complex spatial 
structure of the cloud. The fit of the 2D Ay map is done 
by determining the optical depth along the LoS through a 
3D density structure described by the series. This means 
that parameters describing the spatial structure along the 
line of sight will not be determined by this fit. Instead, a 
library of possible spatial structures is created and used for 
the modelin g. Details of the techniq ue to fit the Ay map are 
described in lSteinacker et all (|2005| ). but here we just sum- 
marize the basic approach and refer readers to that paper 
for details. The dust properties arc fixed for this determina- 
tion to A =0.1 /im silicate spheres, but in later sections we 
re-interpret the obtained optical depths using other grain 
sizes. 

The underlying 3D density structure is described by a 
series of 100 3D Gaussian density clumps with their main 
axes aligned with the Cartesian coordinate axes. The po- 
sition and extent in the plane of sky of each Gaussian (4 
parameters) and their normalization factor (5th parameter) 
are optimized with a, special simulated annealing search tool 
(iThamm et al.|[l994l) . The high number of Gaussians is nec- 
essary because the Ay map is very clumpy with many local 
maxima and structure on small and large scales. Tests with 
lower numbers of Gaussians could not reproduce the main 
features of the map adequately. It must be stressed, how- 
ever, that the original map in Paper I is incomplete (lower 
limit to the extinction toward the northern prestellar core) 
and that it probably suffers from aliasing in the interme- 
diate opacity range (Ay= 20 - 40 mag) because there are 
few background objects. 

In Fig. |U we show the inner 66000 AU of the Ay map 
(left) and the corresponding fitted map using Gaussian 
clumps (right). The underlying Gaussian distribution is 
sketched in Fig.[5]by showing ellipses with the main axes of 
the Gaussians at their corresponding position in the plane 
of sky. The line thickness is a measure of the mass in the 
different clumps. The figure shows that the clumpy struc- 
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Fig. 6. Left: Local optical depth in a cell of the 3D density structure (constant density) as a function of the grain size 
(assuming dust mass conservation). The solid (dashed) line refers to absorption (scattering), and the thick (thin) line 
to a wavelength of 0.55 /xm (3.6 /im). Middle: Optical depth for dust following Eq. JI| as a function of the cell gas 
number density for absorption and scattering at A =3.6 /xm (ao = 6x 1CP 8 m, a = 0.4, nH 2 -,o = 10 10 m~ 3 ). Right: Local 
absorption and scattering optical depth along a typical ray crossing one of the 10,000 3D density data cubes that fit the 
Ay map. 




Fig. 5. Position and main axes of the Gaussian basis func- 
tions used to fit the inner part of the Ay map. The coordi- 
nates x and y are in the plane of sky. 



ture of the Ay map requires a scries of basis functions in 
the inner 10000 AU along the x-axis to reproduce the ba- 
nanashape of the cloud center. In turn, some of the detailed 
spatial structures (e.g. in the lower left part of the Ay map) 
are represented by a low number of basis functions, which 
indicates that the chosen number of 100 is still not enough 
to describe all features of the map. 

The fit of the Ay map provides us with a column den- 
sity structure, but the extent and position of the Gaussians 
along the line of sight are still unknown. To provide a basis 
for the fits of the scattered light images, we have created 



a library of 10000 3D density structures where the spatial 
structure parameters along the LoS were chosen randomly 
given the outer physical borders and the constraint to agree 
with the fit of the Ay map as shown in Fig. 2J 



4. Scattered light image modeling using a dust 
growth model 

4.1. Dust model 

It is expected that in the cold core environment, the parti- 
cles grow to fluffy aggregates, vary in their chemical com- 
position and fo llow an entire size distribution (see, e.g., 
IOssenkopj[l993T) . While there are no direct measurements 
of the dust grain size distribution in dense cores, model 
calculations and laboratory experiments at least can pro- 
vide rough estimates of the mean grain size as a function of 
the g as number density in one free-fall time (|Flower et al.l 
20050. Therefore, we assume that the grains have properties 
approximately described by spherical homogeneous parti- 
cles of one size and chemical composition. The radius of 
the grains a is assumed as a simple powcrlaw of the H2 
number density tih 2 with the spectral index a 




n H2 > n H2 . 
tih 2 < n H . 



(1) 







so that the grains in regions below the threshold density 
n-H 2 ;0 have the size ao. The calculations leading to the best- 
fitting results described in this secti on were performed wit h 
varying values for ao, tih 2 ;0, an d a. lPoppe fc Blum] (1997 ). 
e.g., found a growth index of a = 0.25 (see iFlower et al.l 
I2005L their Fig. 3, plotting the dependence in double- 
logarithmic scale). 

In Sect. 02 we modeled the local optical depth for extinc- 
tion T ex t(X) for A = 0.55 /im toward the observer for each 
cell of our spatial grid. For the radiative transfer model- 
ing of the Spitzer scattered light images, we need to cal- 
culate the optical depth for absorption and scattering at 
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the Spitzer band wavelengths in each cell. If we know the 
grain size and the number density of these grains, n, in each 
cell, we can use existing opacity tables to find the cross sec- 
tions cr ahs / sca (A, a) and calculate the optical depths. To use 
Eq. ([!]), this requires finding the gas density from the de- 
rived local extinction T ext (0.55 pm). 

The optical depth across a Cartesian grid cell with the 
edge length s is 

T ext (0.55pm) = a ext (0.55pm, A) n s, (2) 

if all particles have the size A and no = n(A). Assuming an 
internal density of each homogeneous grain pd, the mass of 
a grain is nid = An A 3 pd/3, and the total dust mass in the 
cell volume is n md s 3 . For a given mass gas-to-dust ratio 
Rg/d, this defines the H2 mass fijj a m# 2 s 3 , hence 

4ttA 3 p d _ T ext (0.55pm) 1 

3mH 2 a ext {U.55pm,A) s 

This allows us to calculate the H2 densities for all cells and 
the corresponding grain size using Eq. ((T|) . The correspond- 
ing number density n{a) = no A 3 a~ 3 results from the con- 
servat ion of dust mass in the cell. Using the iDraine fc Led 
(|l984h opacities, we can then calculate the corresponding 
local MIR optical depth r = a a b s / sca (3.6pm,a) n(a) s 
for the absorption and scattering needed for the radiative 
transfer modeling of the scattered light images. 

The effect of grain growth on optical depth due to ab- 
sorption and scattering is illustrated in Fig. [6] The left 
plot shows the optical depth for absorption and scatter- 
ing through a typical cell of the 3D density structure with 
constant density for A =0.55 /im and 3.6 /im as a function 
of the grain size. The grain numbers have been chosen to 
conserve the dust mass in the cell. For optical wavelengths, 
the change in the ratio of <J a bs and <r sca is less than a factor 
of 4 over a tenfold increase in size. In contrast, at 3.6 /im, 
absorption becomes more efficient by a factor of 3, while 
scattering is amplified by a factor of 300. Even a moderate 
increase of grain size by a factor of 2 will result in about 10 
times more efficient scattering. Correspondingly, the MIR 
wavelength range allows to see larger grains in scattered 
light as scattering becomes efficient again for higher gas 
number densities above about uh 2 -q if the optical depth 
for absorption allows the scattered light to proceed toward 
the observer. The dependence of the optical depth for ab- 
sorption and scattering on the gas number density when 
using Eq. (TT]) is illustrated in the middle figure. In the dens- 
est regions, the optical depth for scattering is larger than 
for absorption at 3.6 (im. The right plot gives the opti- 
cal depths for a typical ray cutting an entire density cube 
consistent with the Ay map of L183 and including grain 
growth. In this figure we see that while we cross clumps of 
high density, scattering becomes stronger, so that scattered 
light should be seen near the densest regions in contrast to 
NIR light, for which more scattering is seen in the outer 
parts. Whether scattered light indeed reaches the observer 
from the inner parts is a complex interplay of background 
radiation, optical depth, and forward scattering, and can 
only be answered by running a radiative transfer code. 

Aside from the change in the scattering efficiency due 
to the dust particle growth, the scattering phase function 
describing how radiation is scattered from one direction to 
another also changes. Figure [7] shows a typical phase func- 
tion in polar coordinates of a homogeneo us spherical sili- 
cate particle at a wavelength of 3.6 pm (jSteinacker et alJ 
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Fig. 7. Phase function in polar coordinates of a homoge- 
neous spherical silicate particle for different particle sizes a 
at a wavelength of 3.6 pm. The H2 number densities cor- 
responding to each a as defined in Eq. |T]) are given in the 
lower right part. 

|2003[ ) placed at the origin. While small ISM particles scat- 
ter almost isotropically at this wavelength, the scattering 
becomes peaked in the forward direction for larger grains. 

This forward beaming can alter or even dominate the 
appearance of the scattering images. For the outer region of 
the core, the interstellar radiation field (ISRF) reaches the 
dust from many directions and is scattered almost isotrop- 
ically. For the central region, mainly forward scattering 
dominates because of grown particles. The interstellar radi- 
ation from the direction of the observer will not be scattered 
in a backward direction. Radiation hitting the core in the 
plane of the sky also has a low efficiency to be scattered by 
90°. Only radiation from the background is scattered effi- 
ciently in a forward direction toward the observer. But this 
radiation has to pass the central parts of the cloud first and 
will be seen only if the optical depth is still so little that 
emerging flux can be detected. For Spitzer, this might actu- 
ally limit the possibilities of observing the growth of grains 
in the regions of highest density. However, with the sen- 
sitivity of the mid-infrared instrument onboard the James 
Webb Space Telescope, the wavelength range > 5 pm. will 
become important for tracing larger grains. 

4.2. Theoretical cloud shine images of L183 and 
dust-modeling results 

The presentation of the image modeling results in this sec- 
tion has two goals: (i) to demonstrate that scattered light 
images based on a density structure in agreement with the 
current extinction data of L183 can reproduce the main 
features of the observed Spitzer images, and (ii) to present 
evidence that only models with grain growth can repro- 
duce the appearance of the central banana-shaped region 
in contrast to models with ISM grains without growth. 

For the first goal, it is sufficient to present one image 
with a distinct set of parameters of the density structure 
model, the dust data, and the external field that shows 
overall agreement. The common approach for the second 
goal would be an automated fit and exploration of the pa- 
rameter space quantifying the agreement on the basis of a 
X 2 -minimization as we used it for the Ay-fit. An automated 
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Fig. 8. Comparison of the three Spitzer images at 3.6, 4.5, and 8 /jm of the inner 66000 AU of L183 (top) with scattered 
light models based on grains growing as a function of density (bottom). The underlying 3D structure of the model images 
is consistent with the measured Ay map. The general pattern of the modeled diffuse emission is similiar while, the flux 
is about a factor of 2 lower in the model. 



fitting procedure varying the many free parameters is pro- 
hibitive, though as radiative transfer calculations applied 
to a clumpy 3D structure like L183 including scattering are 
time-consuming. For the images shown in this section, we 
used 50 3 spatial grid cells and a discretization of the direc- 
tion space in 648 directions leading to a intensity solution 
vector with 8.1 x 10 7 values. An automated fitting proce- 
dure would need to repeat the calculation of this vector 
at least several ten thousand times to optimize the density 
structure. 

Therefore, we present here the results of a nonauto- 
mated exploration based on the 10000 spatial structures 
derived from the Ay map fitting as a representation of pos- 
sible density distributions of the cloud. The quality of the 
fit was decided by eye: reproduction of the central banana- 
shaped region and its maximum in the northern part, the 
radiation depression in the central core, the overall agree- 
ment with the flux level observed with Spitzer, and the 
shape and location of the weak emission patterns in the 
outer region. 

Using the ray-tracing version of our 3D radiative trans- 
fer code (Steinacker et al. 2003), we selected spatial struc- 
tures, which yield images fitting the Spitzer 3.6 /im image 
for varying dust models. Then we further changed both the 
structure and the dust parameters by hand to study the 



influence of each parameter and to find the best-fitting pa- 
rameter set. Finally, we calculated the images at 4.5 and 
8.0 /xm for the same parameter set. 

As an external radiation field, we used the DIRBE 
data provided by the legacy archive microwave back- 
ground data analysis (LAMBDAQ) to take the global 
anisotropy of the ISRF into account. Moreover, a nearby 
star could be another (local) source of anisotropic illu- 
mination. Investigating the stars in the vicinity of LI 83, 
we found one star that could be near or in the cloud. In 
Appendix [Bl we discuss the spectral energy distribution 
(SED), and the resulting stellar properties. A star consis- 
tent with the SED and placed near or even in the outer 
parts of LI 83 would locally influence the scattered light 
image at 3.6 [im up to distances of 2000 AU from the star. 
In all other parts of L183, the ISRF dominates the radiation 
field. Therefore, we left the star out in our consideration, 
like the other stars visible in the IRAC images. 

Figure [5] compares the Spitzer images with the result- 
ing scattered light model images. The upper images show 
the inner 66 000 AU of the cloud at 3.6, 4.5, and 8 /jm as 
observed by Spitzer (see Fig. [I}. Images calculated with- 
out growth of the grains were too low in scattered light 
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flux at 3.6 /im for various ISM grain sizes, and all fluxes 
were below 10 -5 MJy/sr. The lower images show the best- 
fitting model images at 3.6 (left), 4.5 (middle), and 8.0 /im. 
The corresponding parameters of the growth model as de- 
scribed by Eq. |T]) are a=0A for the powerlaw index, a 
grain size in the outer region of ao = 5 x 10~ 8 m, and a 
H2 number density threshold for growth at 3 x 10 10 m -3 . 
In the lower left hand image, the overall emission pattern 
is well- matched with more emission in the northern part. 
The model images show some peaks (e.g. directly below the 
central core), which obviously are not present in the Spitzcr 
images. Investigating the underlying density structure, we 
found that those are attributed to single Gaussian clumps 
and the image peaks could not be removed by changing 
their LoS position or extent. We relate them to either a 
remaining error in fitting the Ay map or to uncertainties 
in the Ay map itself. 

The lower middle 4.5 /im image also shows overall agree- 
ment with the corresponding Spitzcr image. The same flux 
normalization is used as in the 3.6 /im image. The up- 
per part of the central banana-shaped region is visible, 
and the flux ratio in the observed and theoretical images 
F4. 5/j,m/3.6nm is about 0.2. The 3.6 /an Spitzer image shows 
the central core in absorption and also in the 4.5 /urn image 
this dip is weakly seen (compare with Fig. [5]) . In the model, 
the dip is reproduced at both wavelengths. The actual size 
of the dip is not reproduced because of the limited size of 
the computational cells. 

For the 8 /im image, a background flux value was inter- 
polated from regions of low extinction in the Spitzer image 
and assumed as a constant background. Therefore, com- 
paring the image with the Spitzer data has to be done with 
caution: bright regions in the lower right image correspond 
to regions where the background radiation experiences lit- 
tle extinction. The calculations include scattering (which is 
very low even in the presence of larger grains) , but do not 
consider the effect of PAH emission by the cloud itself. 

As absorption dominates at all densities, the image re- 
produces the main features of the Ay map. It also reveals 
where the density structure model was still too coarse: sin- 
gle Gaussian clumps appear where the Spitzer image shows 
more complex absorption structures. The overall shape nev- 
ertheless agrees with the observed data. We also cannot 
exclude the influence of variations in the background. 

The fitting process revealed a prominent ambiguity 
in the obtained best-fitting dust growth parameter sets. 
Comparing the sets, we found that they are based on the 
dominance of larger particles to create scattered light in 
the MIR. As a result, the images were not sensitive to the 
density threshold and the outer particles size, as long as the 
outer particles are smaller than 0.1 /im to avoid substan- 
tial scattered light flux in the outer parts. This ambiguity 
is visible in the growth law defined in Eq. (JTJ) : changing 
the density threshold by a factor / can be compensated 
for by adjusting the size of the outer grains by a factor 
f a to get the same shape of the grain growth. In contrast, 
changing the spectral index had a strong effect on the im- 
ages. Already for a = 0.3, the flux values were too low by 
a factor of 4. 

In general, the flux of the theoretical images is about 
a factor of 2 lower than the observed flux. We show the 
specific flux density of model and observed image in Fig. [5] 
both as flat image and in surface representation to compare 
the actual values. As the stars confused the image substan- 



tially, we have processed the 3.6 /im IRAC image through 
a routine that removes the stars. The routine considers the 
flux distribution of subimages and marks pixels with fluxes 
in the tail of the distribution, then interpolates them from 
the remaining pixels. The emission peak near the bright 
star in the lower left part of the image is an artefact from 
the spread-out PSF of this star. 

While the flux pattern could be reproduced by the 
model, the higher flux could not be achieved by chang- 
ing the grain properties or the spatial structure. Within 
the simple dust model, one might consider that higher flux 
value can be achieved by assuming a stronger growth of 
particles with the bigger grains producing more scattered 
light. The increased extinction, however, damps this radi- 
ation more efficiently so that the emission of these grains 
cannot be observed any more. 

The flux will also change if we consider a more sophis- 
ticated dust model. As an example, if we consider cllipti- 
cally shaped grains with perfect grain alignment to a large- 
scale magnetic field, the phase function, hence the scat- 
tered light flux, will increase by a factor of 2 by assum- 
ing a moderate grain main axes ratio of just 2 (see, e. g., 
IVoshchinnikov fc Farafonovlll993D . 

It can be expected that an automated fitting and/or 
an improved modeling of the Ay map would remove most 
of the obvious peaked artefacts. Further improvement can 
be expected from improving the Ay map in the northern 
banana-shaped region because in the original study (Paper 
I), that part was not mapped with the MIR camera of the 
Infrared Space Observatory (ISOCAM). From mm mea- 
surements, the expected extinction reaches Ay s=s 60 to 70 
mag (being somewhat insecure for technical problems with 
the bolometer). This was beyond the sensitivity of the H-K' 
color-reddening method used to map this part to comple- 
ment the ISOCAM map. Lower limits in the range Ay = 
25-40 mag have been filled in instead, which did not induce 
dense clump fitting in the 3D Ay model. 

In a more realistic model, larger grains might be ex- 
pected to grow as a second population from the original 
seed population. So there is the normal Rayleigh scattering 
by most particles and some more efficient scattering of big- 
ger grains. In the simple model in the paper, all grains are 
moved away from the Rayleigh regime when growth occurs. 

In test runs, we also estimated the limit for the largest 
grains we are able to detect with our method in L183 based 
on our dust model. For 3.6 and 4.5 /im, grains of about 1 
/im in size make the largest contribution to the scattered 
light flux. The flux from larger grains quickly drops because 
of the strong extinction within the dense regions. 

An interesting point is the timescale of the growth pro- 
cess related to the evolution time scale of the core and 
cloud. If the clouds are long-lived, grains might grow enough 
to become a seed population for the subsequent planet for- 
mation process. Th e growth model that we adopted from 
iFlower et al.l (|2005l ) calculated the growth status after one 
free-fall time. We will investigate this question further in a 
forth-coming paper dealing with the growth and maximal 
sizes expected for molecular cloud cores. 

To sum up, we refer to the two goals formulated at the 
beginning of this section: (i) we have demonstrated that 
scattered light images based on a density structure that 
agrees with the current extinction data of L183 can re- 
produce the main features of the observed Spitzcr images. 
Some deviations arise from numerical artefacts or incom- 
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Fig. 9. Comparison of the specific flux density (z-axis) of the 3.6 /zm Spitzer image of the inner 66000 AU of L183 
(right), processed through a star removal routine, with the scattered light model image (left). 



pleteness of the extinction data and are expected to be 
tackled with a better modeling; (ii) the presented images 
use a grain growth model, and indeed the images are very 
sensitive to the threshold density where growth starts, to 
the starting grain size, and to the growth efficiency We find 
that models without dust growth do not produce any mea- 
surable scattered light flux, while a growth from grains with 
a size of about 0.05 /xm starting at an H2 number density 
of 3 x 10 10 m~ 3 can reproduce the emission pattern seen 
by Spitzer in the 3.6 and 4.5 /im filter. 



5. Summary 

In this paper, we have presented Spitzer maps of the dark 
cloud complex L183 which reveal emission at 3.6 and 4.5 
/im and absorption at 8 /im near the densest regions of the 
cloud (Ay > 15). We argued that the emission cannot be 
attributed to thermal black body emission of larger dust 



grains, they are too cold in this environment. The emission 
of stochastically heated particles should originate in the 
layer where most of the interstellar radiation field is ab- 
sorbed. The morphology of this layer is visible in the scat- 
tered light at optical wavelengths and traces the outer con- 
tours of the extinction map, not the inner banana-shaped 
region of L183. Moreover, existing PAH emission models 
show a strong increase in emission from 4.5 to 5.8 /xm, 
which is not seen in the Spitzer images. Finally, a careful 
investigation of the density structure has confirmed that 
an excitation of stochastically heated particles (PAHs and 
small grains) through holes in the clumpy cloud is unlikely, 
because of the optical depth for the interstellar radiation 
at optical and UV wavelengths for almost all investigated 
spatial structure models for L183 (see Appendix A). We 
therefore suggested that this emission is scattered inter- 
stellar radiation. 
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To test this hypothesis, we performed 3D radiative 
transfer calculations producing scattered light images of 
the central dense region where emission is seen. The den- 
sity structures were derived from an automated fit of the 
Ay map for this region (from Paper I) by a series of 3D 
Gaussian number density functions. The missing informa- 
tion along the line-of-sight was dealt with by creating a 
library of possible density structures in agreement with the 
Ay map. 

When using simple silicate dust grains without growth, 
the produced scattered light flux was not sufficient to re- 
produce the observed image. This finding was independent 
of the chosen spatial structure and the dust grain size. The 
model included a realistic direction-dependent ISRF based 
on the D1RBE data, which did not introduce new unknowns 
that could change this view. 

This led us to consider grain models that allow for grain 
growth in the denser parts of the cloud. Assuming a sim- 
ple growth law along the directions of recent laboratory 
measurements, we investigated the changes in the optical 
depth and the scattering direction for photons propagating 
into regions with increasing density containing larger dust 
grains. 

Within this model, we were able to reproduce the main 
features of the Spitzer observation in the 3.6, 4.5, and 8 
/zm bands: (i) the emission arises mainly in the dense re- 
gions, (ii) the emission in the 4.5 /jm band is weaker than 
in the 3.6 /im band and concentrated in the densest upper 
part of the banana-shaped central region, (iii) the central 
region appears in absorption at 8 fim. The flux is smaller 
by a factor of about 2 in the model images compared to the 
3.6 (im flux measured by Spitzer. However, various uncer- 
tainties enter the modeling, namely the uncertainties in the 
Ay map, the errors in the Ay map fitting, the uncertain- 
ties in the extinction properties of the dust, the uncertainty 
in the scattering phase function, the approximations in the 
dust model, and the uncertainty in the ISRF. 

The existence of a new MIR window to study the growth 
of grains in the dense parts of cores is the result of an inter- 
play of several effects: At MIR wavelengths, the moderate 
optical depth allows observing dense regions where growth 
is expected. The MIR scattering efficiency of smaller grains 
is too low to detect them, but larger grains scatter more 
efficiently and become observable again. 

From our results we conclude that this is the reason 
why we see emission in the 3.6 and 4.5 /im Spitzer images 
of LI83 as interstellar radiation scattered by dust particles 
grown in the dense region with respect to common ISM 
particles. Therefore, the Spitzer observations of LI83 are 
the first direct evidence of dust growth in molecular cloud 
cores. In turn, we evaluate the evidence from the extinction 
measurement at MIR and from FIR/mm observations to be 
indirect as there is an ambiguity in explaining the data by 
either assuming larger grains or by modifying the density 
and temperature structure used in the LoS integration. 

In a future paper, we will present a 3D model of the 
dust based on a revised Ay map; and it will be strongly 
constrained by combining all available data from mm to 
NIR. 
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Appendix A: Exploring the applicability of a PAH 
model: Three-dimensional structure variations 
and r-iso-surfaces 

From the Ay map alone, it cannot be excluded that stellar 
emission is approaching the inner parts of LI83 to excite 
PAHs at the locations where emission is seen in Spitzer 
images. Projection effects can hide holes which would allow 
such excitation, and would not to be unlikely given the 
clumpy structure of the core. Becausethe PAHs are excited 
by optical and UV photons, we consider a wavelength of 
0.55 /im as a conservative limit - if this radiation is not able 
to reach the PAHs then radiation at shorter wavelengths 
will not cither. In this analysis, we define "reaching" as 
t < 1 . We performed spatial structure analysis calculations 
to explore the cases for which PAH excitation by the ISRF 
is possible. Varying the extent and position of the Gaussians 
along the linc-of-sight randomly we obtained 10000 spatial 
structures, each of it in agreement with the LI 83 Ay map. 
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Fig. A.l. Minimal optical depth at 0.55 /im r m to the 
central density maximum as a function of the 3D distance 
of the To.55fj, m = 1-surface to the central density maximum. 
The upper and lower plots show spatial structures for which 
the central density maximum is located within and outside 
the inner 30000 AU of the spatial structure, respectively. 
The numbers in the lower right corner indicate the number 
of spatial structures in this category. 



In the following, we would like to examine which of these 
spatial structures would allow for a PAH emission model. 
To characterize each configuration, we calculated the min- 
imal 3D distance of the To.ss^m = 1-surface to the central 
density maximum D m (as an indicator for density holes) 
and the minimal optical depth at 0.55 /im toward the den- 
sity maximum r m (as indicator for the heating of the center 
by the ISRF) . The results are plotted in Fig. IA.ll which re- 
lates both quantities. 

To distinguish the spatial structures with a central den- 
sity maximum close to the front or back parts of the consid- 
ered density cube (which is close to the low-density outer 
cloud part) and those that have the center in the mid- 
dle part, we plotted those spatial structures in the upper 
(lower) panels for which the central density maximum is 
located within (outside) the inner 30000 AU of the spatial 
structure, respectively. In the following, we refer to these 
two cases as central and border maximum spatial structures. 



Their numbers are given in the bottom right corner of each 
plot. 

In the 3.6 fxxa Spitzer image, emission is seen out to 
-Dg=4000 AU measured in the plane of sky from the cen- 
tral density maximum. Assuming that the emission comes 
from the surface of a cylinder with this radius, the ISRF 
has to enter the cloud down to this distance to excite PAHs. 
The dash-dotted line indicates this distance De and the 
1-surface should be closer to the center or ef- 



70.5 



ficient PAH excitation will not be possible. As a conser- 
vative criterium, we assume that minimal distance of the 
To.55/im = 1-surface to the center should be less than De- 

For most border maximum spatial structures and al- 
most all central maximum spatial structures (right of the 
dash-dotted line) this is not the case. Exciting stellar opti- 
cal light cannot reach the inner parts without being heav- 
ily extincted, so that PAH emission can be excluded as an 
explanation for the observed emission. A few border max- 
imum spatial structures reach low r m < 10 values. The 
resulting temperature of the ce ntral density peak would b e 
higher than the 7-9 K found bv lPagani etafl (|2003l . |2004[ ) . 
Thus, spatial structures below the dashed line at r m < 10 
are also not consistent with the observational data. 

We arc therefore left with spatial structures in the up- 
per left hand corner of the two plots in Fig. lA.ll to be can- 
didates for possible PAH emission models. For the upper 
plot, the central density maximum can be located in the 
front part or in the back of the entire spatial structure. If 
it is located in the back, the excited PAH emission has to 
cross the entire spatial structure. Since T^.Q^m is about an 
order of magnitude less than To.s^m, the expected optical 
depth would still be between a few to ten, making the PAH 
emission produced in the front part dominate the emission. 

If it is located in the front, it will be illuminated well, 
hence expected to be seen in emission. Only if there is more 
flux created in the back by the other less dense clumps 
would it be seen in absorption. This region between the 
other clumps and the foreground clump are shielded well 
though, because to excite PAHs, radiation has to cross ei- 
ther the other clumps or the foreground clump. 

Limb brightening in an optically thick clump would pro- 
duce a ring of high surface brightness and a dip in the cen- 
ter. But such an external dense clump should be seen in 
the optical images of LI 83 as a small region with enhanced 
scattered light. But the CFHT I image only shows smooth 
scattered light on scales of the outer r = 1 layer visible in 
the Ay map. Moreover, just a dip is seen but not a ring. 

In view of these arguments, it appears unlikely that the 
central part of LI 83 has a density maximum located near 
the edge of the spatial structure. The remaining few central 
maximum spatial structures in the lower plot of Fig. I A. 1 1 
left from the D m dash-dotted line and with r m > 10 are the 
only spatial structures with holes arranged in a way that a 
PAH emission model would work. 

With 10000 random spatial structures, a parameter 
space of 200 parameters is not well-sampled, but the main 
optical depth effects are dominated by a few ten Gaussians. 
From this analysis, we conclude that a PAH model for the 
3.6 fj,m emission cannot be excluded formally with the avail- 
able data. From our sample of 3D spatial structures con- 
sistent with the Ay map, only a very small fraction would 
allow for PAH emission near the central absorption peak, 
which therefore seems unlikely. 



12 



Steinacker et al.: Direct evidence of dust growth 



Appendix B: A possible star within L183 

Beside the ISRF, another source of illumination could be 
the presence of a nearby star, positioned at ~3' away from 
the reference position. This star, visible in neither the blue 
POSS plates nor the red POSS I plate and hardly visible 
in the red POSS 2 plate, is bright in the infrared (Table 
IB.ip . implying that the star is attenuated by the cloud. We 
started by deriving a set of plausible reddening values (from 
the H and K' difference, we found Ay in the range 10 to 23 
mag). For this range, no normal star (main sequence, giant, 
supergiant) could fit the data. 

Table B.l. Star position and magnitudes 



Position 


R.A. 


Dec 




(J2000) 


15 h 54 m 20.45 s 


-02°54'07.30" 




Band 


Brightness 


Error 


Ref 


(H 


mag 


mag 




R 


(19.7) a 




1 


I 


16.3 


0.1 


2 


J 


10.986 


0.021 


3 


H 


8.915 


0.044 


3 


K 


7.910 


0.016 


3 


3.6 


7.361 


0.002 


4 


4.5 


7.340 


0.002 


4 


5.8 


7.222 


0.003 


4 


8 


7.177 


0.003 


4 


24 


7.4 


0.1 


4 



a USNO-B calibrations are done in the original E band 
plate, which differ s somewhat from the standard R band 
([Monet et alj|2003f ) 

References : (1) USNO-B flMonet et al.ll2003f) (2) u npublished 
data from Paper I, (3) 2MASS (|Cutri et al J l2000h . (4) this 
work. 

We therefore applied a more detailed analysis, by apply- 
ing an automated SEP fitting procedure ba sed on the sim- 
ulated annealing code (|Thamm et al.|[l994l ) that calculates 
the spectrum using an extincted Planck spectrum varying 
the effective stellar surface temperature T, the distance , 
the stellar radius R, and the column density of extinct- 
ing material for ISM Si and C grains. Figure IB. II shows 
the specific spectral energy density \F\ as a function of 
A. The data points are marked as thick crosses, no error 
bars have been considered as the given errors are small (see 
Table [B~Tj) . Comparing with an unextincted spectrum in 
the Ray leigh- Jeans limit (dotted line), no indication of IR 
excess is visible, which would arise from the thermal emis- 
sion of warm dust in a nearby disk and/or shell. 

We tried two ways to fit the spectrum with a stellar 
spectrum that undergoes no extinction. The dashed line 
gives a fit for a Planck spectrum that has the same maxi- 
mum as the observed SED (T» = 2400 K, fl, > 3 x 10 6 km, 
same distance as L183). Cool object spectra fail to repro- 
duce the sharp turnover around the maximum without ex- 
tinction. The dash-dotted lines shows the Planck spectrum 
without extinction that follows the data where the devia- 
tion from the Raylcigh-Jeans spectrum occurs (T* = 7000 
K, i?* > 2 x 10 6 km), but the flux decrease is too weak 
toward shorter wavelengths. Hotter objects cannot repro- 
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Fig. B.l. Specific spectral energy density as a function 
of the wavelength. The observed data points are shown as 
thick crosses. Three fits are displayed without extinction 
are displayed: a spectrum in the Rayleigh-Jeans limit (dot- 
ted), a Planck spectrum having the same maximum as the 
observed SED (dashed), and a Planck spectrum that fol- 
lows the data where the deviation from the Rayleigh-Jeans 
spectrum occurs (dashed-dotted). For the case of extincted 
stellar emission, two spectra are shown: a Planck spectrum 
of a star being extincted by silicate grains (thick solid), and 
by carbon grains (dash-three dotted), respectively. 



duce the sharp turnover either. We therefore concluded that 
extinction is needed to achieve better fits. 

Putting a layer of extincting dust in front, we used the 
SED optimizer to investigate two cases: (i) silicate grains 
(solid thick line) and (ii) carbo n grains (solid dash-three 
dotted line) (jDraine fc Led[l984[ ). The thick solid line gives 
the fitted spectrum for silicate grains, the dash-three dot- 
ted line for carbon grains. Both fits failed to reproduce the 
SED with good accuracy. In these models, the temperature 
was around T» = 7000 K, and the extinction was about 
t =5 to 9 depending on the dust model. It would need a 
more sophisticated dust and/or stellar atmosphere model 
to improve the fits. Given the projected position in the ex- 
tinction map, the optical depth from the Ay map for the 
star when sitting behind the core would be 10 to 15. This 
value is close to the extinction derived from the SED fit- 
ting. We conclude from these data that the star is not in 
front of the cloud. 

We calculated the flux ratio of a star consistent with the 
SED fits at various distances from L183 and the ISRF at 
3.6 /im using DIRBE data. We find that the radiation field 
is dominated by the star up to about a distance of 2000 
AU and then the stellar influence drops quickly due to the 
quadratical decrease with distance. If the star were placed 
near LI 83 or even exactly at the same distance to Earth 
as the density maximum, it still would be located in the 
outer parts of L183 as is evident in Fig. 1, and we would 
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only expect a very local (a few thousand AU) influence on 
the scattered light pattern, while the ISRF dominates for 
the entire banana-shaped region. The potential influence of 
the star being placed near LI 83 is weakened further by its 
having to be near a dense region with grown particles to 
produce scattered light at 3.6 micron. 



